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ABSTRACT 

In the lenticular galaxy NGC 1023 a third population of globular clusters (GCs), 
called faint fuzzies (FFs), was discovered next to the blue and red GC populations 
by Larsen & Brodie. While these FFs have colors comparable to the red population, 
the new population is fainter, larger (-R c ff > 7 pc) and, most importantly, shows clear 
signs of co-rotation with the galactic disk of NGC 1023. We present N-body simulations 
verifying the hypothesis that these disk-associated FFs are related to the young massive 
cluster complexes (CCs) observed by Bastian et. al in M51, who discovered a mass- 
radius relation for these CCs. Our models have an initial configuration based on the 
observations from M51 and are placed on various orbits in a galactic potential derived 
for NGC 1023. All computations end up with a stable object containing 10 to 60% 
of the initial CC mass after an integration time of 5 Gyr. A conversion to visual 
magnitudes demonstrates that the resulting objects cover exactly the observed range 
for FFs. Moreover, the simulated objects show projected half-mass radii between 3.6 and 
13.4 pc, in good agreement with the observed FF sizes. We conclude that objects like 
the young massive CCs in M51 are likely progenitors of the FFs observed in NGC 1023. 

Subject headings: globular clusters: general — galaxies: individual (NGC 1023, M 51) 
- methods: n-body simulations — galaxies: star clusters — galaxies: evolution 



1. INTRODUCTION 



A large number of galaxies show a bimodal color distribution o f their globular cl uster (GC) sys- 
tems, indicating a metal poor blue and a metal rich red GC system (jWest et al.ll2004l . and references 
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there in). lLarsen fe Brodid (|2000l ) studied the GC system of the nearby (9.8 Mpc. ICiardullo et al 
2002) lenticular galaxy NGC 1023 and discovered a clear bimodal color distribution. The average 
effective radius of the blue GCs in NGC 1023 is about 2.0 pc, while the average radiu s of red com- 
pact GCs is slightly smaller (R c s ~ 1.7 pc), as observed also for other galaxies (e.g. lLarsen et al 
200 ll ). The sizes for the blue clusters are as large as -R c ff = 7 pc, while the sizes of the red clusters 
extend up to i? c ff = 18 pc. The objects larger than 7 pc are generally fainter than the red compact 
objects. The spatial distribution of these objects appears to be associated with the lenticular disk 
of NGC 1023. Due to their faint and extended appearance this new star cluster population was 
called "Faint Fuzzies" (FFs). 



Brodie fc Larsenl (|2002l ) performed spectroscopic observations of the GCs in NGC 1023 and 



found that the FFs show a clear sign of co-rotation with the galactic disk, while the compact red GCs 
show no sign of rotation, suggesting that the FFs are indeed a separate population associated with 
the disk of NGC 1023. They averaged all FF spectra to achieve a high e nough signal-to-noise ratio to 
estimate the age of these objects that appears to be older than 7-8 Gyr. iBurkert et al.1 (j2005l ) further 
analyzed the distribution and radial velocities of the FFs in NGC 1023 and suggested that they were 
the remnants of a past gravitational interaction formin g a ring-like structure. Further observations 
also discovered FF s in the lenticula r galaxy NGC 3384 (jBrodie &: Larsenll2002l ). the dwarf irregula r 
galaxy NGC 5195 (|Lee et al.ll2005l ). a number of galaxies of t he Virgo cluster (|Peng et al.l 12000 ) . 
and the lenticular galaxy NGC 1380 (|chies-Santos et al.|[2007l ). 



Observations have shown that young star clusters a r e ofte n not isolated, but are part of larger 
structures called cluster complexes (CCs). iBastian et al.l (j2005l ) observed massive young star cluster 
complexes in the disk of the spiral galaxy M51 and found a scaling relation between the masses and 
the radii of the CCs that is comparable to the relation of the progenitor giant molecular clouds. 
This relation constrains the parameter space for the CC models in our simulations. 



Kroupal (|1998l ) studied t he dynamical evolution of C Cs observed to be forming in the inter- 



acting Antennae galaxies, and iFellhauer &: Kroupa I (|2002l ) performed the first N-body simulations 
indicating that FFs can be formed by merging of star clusters in CCs. The work presented in this 
paper is a follow-up study of these early computations. The simulations were perfor med to test the 



hypot hesis that the FFs are the remnants of the young massive CCs as observed by IBastian et al 
d2005l ) in M51. 



In Section [21 we describe the method and the parameters used for the calculations. In Section 
El we present the results of the simulations, which will be discussed in Section HI 
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NUMERICAL SIMULATIONS 
2.1. Method 



The simulations were performed with the particle-mesh code SUPERBOX (jFellhauer et al 



200(3 ) . The code solves the Poisson equation on a system of Cartesian grids. The code has a 
hierarchical grid structure. The local universe is covered by a fixed coarse grid which contains the 
orbit of the CC around the center of NGC 1023. The CC orbits in an analytical galactic potential 
(disk+bulge+halo). In order to get good resolution of the star clusters two grids with high and 
medium resolution are focused on each star cluster following their trajectories. The individual 
high resolution grids cover an entire star cluster, whereas the medium resolution grids of every 
star cluster embed the whole initial CC. All grids contain 64 3 grid cells. The coordinate system is 
chosen such that the disk of NGC 1023 lies in the x-y-plane. 



2.2. Basic Parameters 



2.2.1. An Estimate of the Gravitational Potential of NGC 1023 



NGC 1023 is a disk-dominated SBO-galaxy at a distance of 9.8 Mpc (jCiardullo et al.l 12002 ) in 
the constellation Perseus. It shows a well behaved rotation in th e inner part, while the ob served 
radial velocities in the outer regions appear to be more random (INoordermeer et al.ll2008l ). This 
behavior is most l ikely related to a past in teraction, which transformed the galaxy into an SB0 
lenticular galaxy. INoordermeer et al.1 (120081 ) argue that the photometry of the disk is consistent 
with a passive fading since 1 Gyr after the galaxy lost its gas. 

However the exact rotation curve is not required for the simulations because a good guess of the 
potential is sufficient to study the evolution of CCs 5 to 8 Gyr ago. In our simulations NGC 1023 is 
represented by an analytical potential, which consists of a disk, a bulge, and a halo component. For 
objects moving on circular orbits within the galactic plane, the results of the simulations depend 
mainly on the total force at a given distance, but very weakly on the exact choic e of bulge, disk, and 
halo p arameters. We use the observed velocities in the center of the galaxy from lSimien Prugniel 
(|1997l ) to constrain the rotation curve in the region of i? ga i < 4 kpc and assume a flat rotation 
curve for the outer part (Fig. [[]). 



The disk is modeled by a Miyamoto-Nagai potential ([Miyamoto &: Nagailll975l ) 



®disk(R, z) 



GM d 



(1) 



R 2 + (a d + Jzi + b 2 d r 



with M d = 5.5 • 10 iU M , o d = 2.8 kpc, and b d = 0.2 kpc. 
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The bulge is represented by a Hernquist potential (jHernquistJ 1 199C1 ) 

GM b 



bulge 



r + a b ' 



with M b = 0.7 • 10 1U M , and a h = 1.1 kpc. The halo is a logarithmic potential 

$halo(r) = \ v l Mr 2 + r halo)i 
with vq = 170 km s , and rhaio = 5.0 kpc. 



(2) 



(3) 



For this project, we neglected heating sources from non-axisymmetric components like Giant 
Molecular Clouds (GMC) and spiral arms. Th is is a good approximation for compact objects with 



masses of about 10 5 M Q . iGieles et al.l ( 20061 ) demonstrated that star clusters with such masses 



cannot be destroyed by individual encounters with GMCs, if the extended nature of GMCs is 
taken into account. Close encounters between GMCs and massive star clusters act as tidal heating 
on the clusters. The corresponding disru ption time is much larger than the integration time of 
the simulations presented in this paper. IGieles et al.l (120071 ) analyzed the effect of spiral arms 
passages on the evolution of star clusters and found for massive clusters as considered in this paper 
a disruption time of the order 100 Gyr. 



2. 2. 2. Initial Configuration of the Cluster Complexes 



Bastian et al.l (|2005l ) derived a mass-radius relation for the CCs in M51 that constrains the 
parameter space for our CC models. They found 11 CCs with ages younger than 10 Myr in the 
disk of M51 that have sizes between 85 — 240 pc, and cover a mass range of 0.3 - 3 • 10 5 M Q . Most 
of these CCs show a massive concentration of clusters potenti ally merging in thei r centers and a 
couple of isolated clusters in their vicinity (see fig. 2 and 3 in iBastian et al.l 120051 ) . The size of a 
CC is defined as the point where the measured color is equal to that of the background stars in 
the disk. This radius corresponds to the cutoff radius of our models. The spatial distribution of 
the CCs is connected to the spiral arms of the disk. Most complexes are found at the outer edges 
of the spiral arms. The CCs show a mass-radius relation comparable to the one of the progenitor 
giant molecular clouds. 

We choose the parameters of CCs in the simulations such that the models cover the two 
extremes (CC_H and CC_L) and a point in-between (CC_M) of the Bastian relation (Fig. [2]). 



The star clusters building up the CCs in our simulations are Plummer spheres (|Plummerlll911 
Kroupal l2008l ). The density distribution of the spheres is truncated at i? cu t = 5i? p i. Each star 



cluster consists of Nq C = 100 000 particles. The CC models are diced according to a Plummer 
distribution with ma sses and radii t aken from the Bastian relation (Fig. [2]). Our models resemble 
the CCs observed by lBastian et al.l (|2005l ) with a high density concentration in their centers. 
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Bastian et al.l (|2005l ) estimated a ver y high star formatio n efficiency within the CCs, using th e 



CO-to-H2 conversion factor for M51 from lBoselli et al.1 ((2002) and a Salpeter IMF (|Salpeterlll955l ). 
Due to the high st ar formation efficiency, gas expulsion has only little influence on the evolution of 
the merging CCs (jFellhauer fc Kroupa II2005I ). Consequently, the distorting effect of gas expulsion 
was not considered in the simulations. 



3. RESULTS 

We have carried out 25 different numerical simulations to study the influence of varying initial 
CC conditions and orbital parameters. The computations are presented in the following way. In 
Section [3. II we investigate the evolution of CCs in three mass regimes (CC_H, CC_M, CC_L) of the 
Bastian relation. In Section [JO] we consider a maximal model with 80 star clusters. For comparison, 
computations without an external gravitational field are presented in Section 13.31 



3.1. Configurations Covering the Cluster Complex Mass Range 

We investigate the evolution of a CC configuration in three different mass regimes (Fig. [2]) for 
circular orbits at galactic distances of i? ga i = 2, 3, 5, 8, 12 kpc (Section 13.1.1)1 and an eccentric 
orbit (Section I3.1.2p . The number of star clusters comprising the CC is fixed to Nq C = 20 in all 
models. In order to get comparable starting conditions for the models CC_H, CC_M, and CC_L we 
introduce a parameter 

R SC 

a = -2- (4) 
" f>CC ' v 1 

where and are the Plummer rad ius of a single star cluster and the Plummer radius of 



the CC, respectively (jFellhauer et al.ll2002l ). The parameter a is a measure of how densely the CC 
is filled with star clusters for an equal number Nq C of star clusters. In general high values of a 
accelerate the merging process because the star clusters already overlap in the center of the CC, 
whereas low values hamper the merging process. We choose a medium value of a = 0.08, which 
leads to a reasonable range of the Plummer radii of the individual star clusters. The grids were 
also scaled to get the same relative resolution for all models. The chosen initial parameters are 
displayed in Table [TJ 



3.1.1. Simulations on Circular Orbits 

We have performed 15 simulations on circular orbits holding the initial star cluster configura- 
tion and the number of star clusters in the complex fixed. The CC masses were restricted to the 
high- and low-mass end (CC_H_20 and CC_L_20) of the Bastian relation and a value in-between 
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(CC_M_20). The orbital parameters are given in Table [2j The circular velocities, v C i TC , at the dif- 
ferent galactic distances, -R ga i, were derived from the rotation curve of NGC 1023 (Fig. [1]). The 
orbital periods lie between 70 and 350 Myr. 

The timescale of the merging process is very short. The merger object forms within a few 
crossing times of the CC. In all cases the merging leads to a stable final object. The properties 
of the merger objects at t = 5 Gyr are displayed in Table [3l The number of merged star clusters 
varies from 10 to 16. The number increases with galactic distance as the influence of the tidal field 
becomes weaker. 



In all cases we calculate surface density profiles and fit these by a King profile (jKingill966l ). 
Three exemplary final merger objects after 5 Gyr are shown in Fig. [3] as contour plots on the 
x-y-plane and as corresponding surface density profiles. 

The tidal radii vary from Rt = 18.8 pc for simulation CC_L_2_20 at R ga _\ = 2 kpc to Rt = 88.9 pc 
for CC_H_12_20 at R gSj \ = 12 kpc. The enclosed mass of the merger objects is defined as the mass 
within the tidal radius. After 5 Gyr 20 to 60 % of the initial CC mass is bound to the merger object. 
The half-mass radius is the radius of the sphere, where half of the mass is enclosed. However, as 
observers can only derive a projected half- mass radius, we calculate also a projected value defined 
as the projected radius in the surface density plot within which half of the mass is included. The 
projected half-mass radius is slightly smaller than the three-dimensional half-mass radius (Table [3]) 
and corresponds to the observed effective radius, -R c ff • 

The enclosed masses and half-mass radii are increasing with galactic distance R ga \, leading to 
small compact merger objects at R ga \ = 2 kpc and more extended massive ones at i? ga i = 12 kpc. 
Comparing the different models CC_H_20, CC_M_20, and CC_L_20 the tidal field has the largest 
impact on model CC_H_20 impeding the merging process. An estimate of the influence of the tidal 
field is given by the parameter 

■nCC 

P= Rp (5) 



(jFellhauer et al.ll2002l ). which is the fraction of the cutoff radius R^ of the CC and its tidal radius 



rCC rp^g i ar g er tli e value of (3, the larger th e impact. The tidal radius is estimated by using the 



Jacobi limit for a satellite on a circular orbit (jBinnev &: Tremaindll987l ). 

(6) 



p>CC _ I Mac y p 



i 



3-v'i (R gal ) J n gaV 

The values for the initial tidal radii of the CCs and the corresponding parameter (5 of the individual 
models are shown in Table [2j The /3-values are decreasing from model CC_H_20 to CC_L_20 and 
from small galactic distances of R g& \ = 2 kpc to large distances of i? ga i = 12 kpc. Therefore model 
CC_H_20 suffers a higher percental mass-loss than simulations CC_M_20 and CC_L_20 and for each 
simulation the impact of the tidal force on the CC is larger for small galactic distances i? ga i- 
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3.1.2. Cluster Complexes on an Eccentric Orbit 

We performed 3 simulations on eccentric orbits to test the influence of the orbit on the prop- 
erties of the final merger object. The orbit was chosen such that the distance of the CC to the 
center of NGC 1023 varies between 3 kpc and 8 kpc. The orbital period is about 100 Myr. Again 
the computations were performed for the three different mass regimes (Fig. [2]). Table [3] shows the 
properties of the final merger objects. 

For simulation CC_H_ecc_20 the value for the enclosed mass lies between the values of the cir- 
cular orbits at i? ga i = 2 kpc and R g& \ = 3 kpc. The half-mass and tidal radii agree more with those 
obtained for the 2 kpc circular orbit. The properties of the final object of simulation CC_M_ecc_20 
and CC_L_ecc_20 are similar to those of the circular 3 kpc orbit. Thus, for the eccentric orbit the 
peri-galactic distance determines the overall properties of the final merger object. 



3.1.3. Different Initial Configurations 

For comparison we also performed simulations for two additional initial configurations in the 
medium mass regime (CC_M_5_c2_20 and CC_M_5_c3_20). Different seeds were used for the random 
number generator to generate two different initial configurations having the same input parameters 
as CC_M_20 (Table [I]). The initial configurations vary considerable due to the low number of 
clusters constituting the CC. The models were placed at a distance of i2 ga j = 5 kpc on a circular 
orbit. 

The two configurations also lead to a fast merger within the first few crossing times. After 5 Gyr 
the final objects have masses of 2.5- 10 4 M (CC_M_5_c2_20) and 2.9- 10 4 M (CC_M_5_c3_20). The 
corresponding projected half-mass radii are 5.5 pc and 5.6 pc (TableEJ. In the case of configuration 
CCJVL5_c2_20 16 star clusters have merged. Whereas for configuration CC_M_5_c3_20 the number 
of merged star clusters is only 13 although its final mass is larger. The final mass of the merger 
object does not scale with the number of merged star clusters as the newly formed object looses 
considerable mass on its orbit around the center of NGC 1023. The final parameters of the merger 
objects are slightly smaller than those of the corresponding model CC_M_5_20 (see Sect. I3.1.ip . 
While the enclosed mass and the half-mass radius of a merger object depend considerably on the 
initial configuration, the general order of magnitude stays the same. 

3.2. A Model With a Large Number of Star Clusters 

In this section we consider a maximal CC with Nq C = 80 star clusters. The parameters for 
the computation are presented in Table [TJ The general parameters for the CCs are taken from 
model CC_H_20. Due to the larger number of star clusters, the individual star clusters have lower 
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masses, which are comparable to the star cluster masses of simulation CC_M_20. The complex is 
placed only on one circular orbit at a distance of R ga _\ = 5 kpc, because this computation is very 
time-consuming. 

The initial star cluster distribution is displayed in Fig. H] as a contour plot on the x-y-plane. 
The distribution is quite extended and has a concentration of star clusters in its center as observed 
for CCs in M51 by Bastian et al. (2005). Within the first 100 Myr a massive merger object forms 
and the unmerged star clusters align along the orbit (Fig. [5]). There is still a close companion star 
cluster as a satellite of the merger object. The time evolution of the merger object is shown in Fig. 
[6] as contour plots. The snapshots were taken at t = 0.5, 1, 3 and 5 Gyr. At t = 500 Myr (top 
panel) the merger object has already formed. It is very extended and shows prominent tidal arms. 
The major mass loss occurs within the first Gyr. Later the merger object becomes more compact 
and its tidal arms are less pronounced. After 5 Gyr the merger object has reached an almost stable 
state. The mass loss is small and the tidal radius hardly changes any more. 

Table [3] lists the properties of the final merger object which resulted from the merger of 32 star 
clusters. Its mass is only 9 % of the initial CC mass and its projected half-mass radius is 8.5 pc. 

3.3. Cluster Complexes without External Gravitational Field 

For comparison we performed simulations without an external gravitational potential (i.e. 
R t = oo) for the three models CC_H_20, CC_M_20, CC_L_20, and the maximal configuration with 
80 star clusters, CC_H_80. In case of the three computations with varying CC masses 17 out of 20 
star clusters merge within the first 500 Myr. After 5 Gyr nearly all star clusters have merged and 
formed an object with projected half-mass radii of 22.5 pc (CC_H_inf_20), 12.7 (CC_M_inf_20) and 
9.0 (CC_L_inf_20) which contain more than 90 % of the initial CC mass (Table [3]). 

For the simulation with 80 star clusters (CC_H_inf_80) the merging process is considerably 
slower. In the first 500 Myr 46 star clusters merge. The newly formed merger object is surrounded 
by several unmerged star clusters which are found preferentially in its close vicinity. During the 
next 500 Myr another 23 star clusters merge into the merger object. After 5 Gyr 79 star clusters 
have merged resulting in a final object with a projected half-mass radius of 45.6 pc and a mass of 
3.9 ■ 10 5 M corresponding to 97.5 % of the initial CC mass. 

The individual star clusters in CC_H_80 have basically the same relative velocities as in 
CC_H_20, but four times lower masses. This reduces the probability of merging considerably. 
The resulting longer merging timescale provides the gravitational field more time to tear clusters 
off the complex. This explains the lower number of merged star clusters and the relatively low mass 
in the case of the corresponding model within a tidal field. 
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4. DISCUSSION AND CONCLUSION 

We summarize the major properties for all computations in Fig. [7J The models CC_H_20, 
CC_M_20 and CC_L_20 were placed on circular orbits at 5 different distances whereas the big 
configuration CC_H_5_80 with 80 star clusters and the two additional configurations were only run 
on a circular orbit at i? ga i = 5 kpc. The 3 computations on an eccentric orbit are indicated by 
dashed lines ranging from 3 to 8 kpc. We plotted the enclosed mass and the projected half-mass 
radius, R e s, against the galactic distance from the center of NGC 1023. 

The enclosed masses were converted into ^-magnitudes to allow for direct comparison with 
the observed data, using the formula 



V = Mi 



V, solar 



2.5 • log(Af FF /x) + 29.97, 



(7) 



where Mv, S oiar = 4.83 mag is the absolute solar V magnitude, Mff the mass of the FF, x the 
mass-to-light ratio a nd the value of 29.97 the distance modulus derived from the planetary nebula 
luminosity func t ion (Ciardullo et al.ll2002l ). We use a mass-to-light ratio of x — 2.3. 9jS observed by 
Prvor Mvlanl (jl993l ) for 56 Globular Clusters in the Milky Way. No observational constraints on 
the M/L-ratio of FFs are available. 

Fig. EH shows the enclosed mass as well as the corresponding V-magnitudes versus the galactic 
distance R^]. The V-m agnitudes of the merger objects cover a range from 22.4 mag to 25.6 mag. 
Larsen fe Brodie (200C) detected FFs in the luminosity range 21.4 mag to 25 ma]], with a median 
of 23.6 mag, but only one FF is brighter than 22.5 mag. The luminosities of the simulated merger 
objects are therefore in very good agreement with the observed luminosities of FFs. 

Figure Wp shows the effective radius of the merged CCs versus galactic distance. The effective 
radii for all simulation s range between 3.6 and 13.4 pc. The effective radius criterion of R c s > 7 pc 
from lLarsen Brodid (12000) is added as a horizontal line. The merger objects of model CC_H_20 
(circular and eccentric) and the big simulation CC_H_5_80 with 80 star clusters lie clearly above the 
radius limit whereas the computations of model CC_M_20 exceed the limit only at galactic distances 
of -Rgai > 8 kpc. The lowest mass merger objects have effective radii below 7 pc. However, the 
effective radius criterion of R c g > 7 pc is not to be taken as a definite limit but more as a reference 
value, because below R^ = 7 pc it is no t possi ble to distinguish FFs from compact red GCs. 80 % 
of the FFs observed bv lLarsen fc Brodid (|2000i ) have effective radii below 14 pc. The median of the 
observed effective radii is 10.7 pc. Our results are therefore in good agreement with observations. 

The simulations presented in this paper demonstrate that CCs from the high-mass end of 
the Bastian relation evolve into FF-like objects. Our formation scenario suggests that FFs are 
the remnants of merged CCs formed in spiral arms. NGC 1023 has a prominent bar that extents 



Larsen fc Brodid (|200oh did not consider objects below 25 mag. Their approximate 50 % completeness limit 



24.5 mag. 
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up to about 3 kpc (IDebattista et al.l 12002). In barre d spiral galaxies, spiral arms usually start 
where the bar ends (e.g. lElmegreen fc Elmegreenlll982l ). The absence of FFs in the inner 3 kpc of 
NGC 1023 is therefore consistent with the formation scenario presented in this paper. The Milky 
Way, which is also a barred galaxy with a comparable size to NGC 1023, has most of its giant 
molecular cl ouds in the disk between Galactic radii 3 to 8 kpc, with a peak between 4 and 5 kpc 
(|Damelll993l ). The mass-radius relation found for young massive CCs in M51 is comparable to the 
corresponding relation for giant molecular clouds in M51 (jBastian et al.ll2005l ). Consequently, the 
radial distribution of FFs should in general follow the radial distribution of giant molecular clouds, 
but modified by the results of Fig. i.e. lower masses and smaller effective radii for decreasing 
galactic distances. In NGC 1023 FFs were found between projected galactic radii 3 kpc < R ga \ < 
10 kpc, with a clear peak at a projected radius of 4.5 kpc. The observed radial distribution of FFs 
is therefore in good agreement with the expected distribution from our model. 

In summ a ry, th e simulated merger objects based on the observed mass-radius relation from 
Bastian et al.1 (|2005l ) resemble the observed parameters mass/luminosity, size and spatial distri- 
bution of FFs in NGC 1023 very well. The merging of young massive CCs is therefore a likely 
formation scenario for FFs. 



The measured velocities from Brodie La 



FFs with the disk of NGC 1023. Burkert et al 



sen! ( 12002 ) confirm the general co-rotation of the 



20051 ) argued that the annular distribution of FFs 



on the sky and some deviations from the rotation curve of NGC 1023 indicate that the FFs form 
a ring-like structure with a mean radius of 5 kpc potentially formed during a past tidal encounter. 
In their model, FFs are on highly eccentric orbits spanning the entire range between 3 and 8 kpc. 
In our model, FFs are associated with the former spiral arms of NGC 1023, as the CCs in M51 are 
associated with the spiral arms of the disk. FFs should therefore follow the rotation of the disk of 
NGC 1023, which was probably disturbed during a past interaction that transformed NGC 1023 into 
a lenticular galaxy. However, du e to the very l o w lum inosities, the observationa l uncertainties of 
the ve locities are relatively large. iBurkert et al.l (|2005l ) re-observed two FFs from lBrodie fc Larsen 
( 20021 ). The new measurements of these clusters yield 539 ±21 and 676 + 13 kms -1 , while the 
previous measurements found 514 ±8 and 725 ± 17 kms -1 , respectively. These differences are too 
large for definite statements about deviations from circular orbits. 



Brodie Larsenl (|2002l ) averaged all FF spectra shifted to zero velocity to achieve a high 



enough signal-to-noise ratio to estimate the age of these objects that appears to be older than 7-8 
Gyr. They used the same integrated spectrum to derive a mean metellicity of [Fe/H] = -0.58+0.24 
and a mean alpha-to-iron ratio [a/Fe] between +0.3 and +0.6 compared to solar values. Enhanced 
alpha-to-iron ratios are usually associated w ith rapid, burs t -like s tar formation on short time scales. 
Indeed, the 5 to 8 Myr old CCs observed by lBastian et al.1 (120051 ) show a low age spread, indicating 
that the stellar mass of the order of 10 5 M Q has been formed in a burst-like event. 

The FFs need to be further analyzed with respect to additional and more precise kinematical 
data to allow for a reliable estimate of their orbits. In addition, age estimates, metallicities, and 
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alpha-to-iron ratios are needed for individual FFs. The scatter and the variation of these parameters 
with galactic distance will shed light on the FF evolution and further test whether the discussed 
formation scenario is valid. 

The work of this paper was supported by DFG Grant KR 1635/14-1. We thank S0ren Larsen 
for helpful discussions and access to their FF-data. 
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R gal [ k P C ] 

Fig. 1. — Rotation curve of a model of NGC 1023 used for the numerical simulations: the circu- 
lar velocity is derived from the chosen analytical potential consisting of a Miyamoto-Nagai disk 
potential (dashed line), a spherical bulge potential (dashed dotted line) and a logarithmic halo 
(solid line). Summing up the three components yields the total potential. The diamond symbols 
represent the data available for the inner part (Simien & Prugniel 1997). 
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Fig. 2. — The mass versus radius relation for CCs in M51 (Bastian et al. 2005). The 11 complexes, 
indicated by diamonds, have sizes between ~ 85 and ~ 240 pc and cover a mass range of 0.3 — 3- 10 5 
M Q . The size of the complex is defined as the radius where the CC is still clearly distinguishable 
from the background light of the galaxy. It is the cutoff radius for our models, which cover the 
high- and low-mass end of the Bastian relation (CC_H and CC_L) and a point in-between (CC_M). 



-15- 




Fig. 3. — Left: Contour plots at t = 5 Gyr on the x-y-plane for the computations CC_H_5_20, 
CC_M_5_20 and CC_L_5_20 at a galactic distance of i? ga i = 5 kpc. The lowest contour level 
corresponds to 5 particles per pixel. The pixel size is 4 pc 2 . This yields 0.25 M© pc" 2 (CC_H), 
0.0625 M pc" 2 (CC_M) and 0.025 M pc" 2 (CC_L). The contour levels increase further by a factor 
of 3. Right: Surface density profile corresponding to the contour plots of CC_H_5_20, CC_M_5_20 
and CC_L_5_20. The profiles are fitted by King models. The dashed vertical lines denote the tidal 
radii. 
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Fig. 4. — Contour plot on the x-y-plane of the initial distribution of the big computation CC_H_5_80 
with 80 star clusters. The lowest contour level corresponds to 5 particles per pixel. The pixel size 
is 4 pc 2 . The other contours increase by a factor of three. 
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Fig. 5. — Illustration of the spatial distribution of the star clusters of computation CC_H_5_80 after 
500 Myr. The CC moves on a circular orbit around the center of NGC 1023 (cross). It consists 
of a large merger object and a leading and trailing arm of about 50 unmerged star clusters. The 
circles indicate the positions but not the sizes of the star clusters. 
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Fig. 6. — Time evolution of the merger object in simulation CC_H_5_80 at a galactic distance of 
Rgui = 5 kpc. Left: Contour plots on the x-y-plane displayed at t = 0.5, 1, 3, 5 Gyr. The lowest 
contour level corresponds to 5 particles per pixel. The pixel size is 4 pc 2 which corresponds to 
0.0625 M pc~ 2 . The contour levels increase further by a factor of 3. Right: Surface density 
profiles corresponding to the contour plots. The profiles are fitted by King models. The dashed 
vertical lines denote the tidal radii. 
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Fig. 7. — a: Enclosed mass M enc i versus galactic distance i? ga i for all computations. Simulation 
CC_H_20, CC_M_20 and CC_L_20 were performed for different distances of R gsi = 2, 3, 5, 8 and 12 
kpc whereas the big computation CC_H_5_80 with 80 clusters and the two additional configurations 
of CC_H_5_20 were only run at a distance of i? ga i = 5 kpc. The eccentric orbit models cover a radius 
range of 3 - 8 kpc and are therefore plotted as dashed lines, b: Effective radius R e d against galactic 
distance R ga _\ for all simulations a s in the left diagram. T he horizontal line at 7 pc indicates the 
lower size limit for FFs defined bv lLarsen Brodid (|2000l ) . 
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Table 1. Initial Cluster Complex and Star Cluster Parameters 



Cluster Complex Star Cluster 



Model 


wCCa 


M cc b 
(M ) 


pCCc 
^cut 

(pc) 


-"pi 
(PC) 




M sc f 
(M ) 


K-SCg 
■"-cut 

(PC) 


(pc) 


CC_H_20 


20 


4 • 10 5 


250 


50 


100000 


2 • 10 4 


20.0 


4.0 


CC_H_80 


80 


4 • 10 5 


250 


50 


100 000 


5 • 10 3 


10.0 


2.0 


CC_M_20 


20 


1 • 10 5 


141.5 


28.3 


100 000 


5 • 10 3 


11.5 


2.3 


CC_L_20 


20 


4- 10 4 


95.5 


19.1 


100 000 


2 • 10 3 


7.5 


1.5 



a Number of star clusters comprising the CC. 

b Initial CC mass. 

c Cutoff radius of the CC. 

d Plummer radius of the CC. 

c Number of star cluster particles. 

f Initial mass of star cluster. 

g Cutoff radius of star cluster. 

h Plummer radius of star cluster. 
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Table 2. Orbital Parameters 









CC_H_20 


CCJVL20 


CC_L_20 


-Kgal 


?; ■ b 

"circ 


rp C 

orb 


i? CCd 











(kpc) 


(km s _1 ) 


(Myr) 


(PC) 




(pc) 




(pc) 


2 


174 


70.4 


42.1 


5.9 


26.5 


5.3 


19.6 4.9 


3 


202 


91.1 


50.0 


5.0 


31.5 


4.5 


23.2 4.1 


5 


220 


139.7 


66.5 


3.8 


41.9 


3.4 


30.9 3.1 


8 


219 


224.3 


91.3 


2.7 


57.5 


2.5 


42.4 2.3 


12 


212 


348.5 


122.4 


2.0 


77.1 


1.8 


56.8 1.7 



a Distance to the center of galaxy NGC 1023. 

b Circular velocity obtained from the rotation curve (Fig. [T]). 

c Orbital period. 

d Tidal radius of CC at t = (Eq. 6). 
e /3-parameter at t = (Eq. 5). 
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Table 3. Results of the 25 Computations after 5 Gyr 



Model 


E i a 

1 l gal 


TV™ b 


M i c 


M i 


V d 




p f 

-n-eff 


1 H 










(%) 




(vc) 






CC_H_2_20 


2 


10 


77923.6 


19.5 


23.5 


8.6 


7.2 


37.6 


CC_H_3_20 


3 


12 


94180.2 


23.5 


23.3 


9.5 


7.8 


41.7 


CC_H_5_20 


5 


12 


116021.0 


29.0 


23.0 


12.3 


10.1 


52.0 


CC_H_8_20 


8 


14 


178244.0 


44.6 


22.6 


15.3 


12.4 


64.3 


CC_H_12_20 


12 


15 


207737.0 


51.9 


22.4 


17.1 


13.4 


88.9 


CC_H_ecc_20 


3-8 


12 


85120.6 


21.3 


23.4 


9.3 


7.4 


38.2 


CC_H_inf_20 


oo 


19 


362804.0 


90.7 


21.8 


29.6 


22.5 


_ 


CC_H_5_80 


5 


32 


35565.1 


8.9 


24.3 


10.6 


8.5 


39.2 


CC_H_inf_80 


oo 


79 


390052.0 


97.5 


21.7 


60.8 


45.6 


_ 


CC_M_2_20 


2 


12 


24366.5 


24.4 


24.7 


5.7 


4.7 


24.4 


CC_M_3_20 


3 


12 


23297.8 


23.3 


24.8 


6.0 


4.9 


26.4 


CC_M_5_20 


5 


13 


34793.1 


34.8 


24.4 


8.3 


6.8 


35.3 


CC_M_8_20 


8 


15 


50045.4 


50.0 


24.0 


10.0 


8.3 


43.5 


CC_M_12_20 


12 


16 


58413.3 


58.4 


23.8 


10.5 


8.3 


58.6 


CC_M_ecc_20 


3-8 


12 


23198.6 


23.2 


24.8 


6.2 


4.9 


26.1 


CC_M_5_c2_20 


5 


16 


24507.4 


24.5 


24.7 


6.9 


5.5 


29.7 


CC_M_5_c3_20 


5 


13 


28668.2 


28.7 


24.6 


6.9 


5.6 


30.4 


CC_M_inf_20 


oo 


19 


90740.5 


90.7 


23.3 


16.6 


12.7 




CC_L_2_20 


2 


13 


11360.3 


28.4 


25.6 


4.4 


3.6 


18.8 


CC_L_3_20 


3 


13 


11879.8 


29.7 


25.5 


4.9 


4.1 


20.9 


CC_L_5_20 


5 


14 


16422.8 


41.1 


25.2 


6.3 


5.1 


25.8 


CC_L_8_20 


8 


14 


20470.2 


51.2 


24.9 


6.5 


5.1 


30.3 


CC_L_12_20 


12 


16 


24679.0 


61.7 


24.7 


7.7 


6.1 


43.9 


CC_L_ecc_20 


3-8 


13 


12491.0 


31.2 


25.5 


4.8 


3.9 


21.1 


CC_L_inf_20 


oo 


20 


37643.0 


94.1 


24.3 


11.7 


9.0 





a Distance to the center of galaxy NGC 1023. 
b Number of merged star clusters. 
c Enclosed mass of merger object. 
d V-magnitudes (Eq. 7). 
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c Half-mass radius of merger object. 

f Effective radius, i.e. the projected half-mass radius of the merger object. 
g Tidal radius of merger object obtained from King fit. 



